Rare White dwarf stars with carbon atmospheres 
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ABSTRACT 

White dwarfs represent the endpoint of stellar evolution for stars with initial 
masses between approximately 0.07 M and 8-10 M Q , where M Q is the mass of 
the Sun (more massive stars end their life as either black holes or neutron stars). 
The theory of stellar evolution predicts that the majority of white dwarfs have 
a core made of carbon and oxygen, which itself is surrounded by a helium layer 
and, for ~80 per cent of known white dwarfs, by an additional hydrogen layer 1-3 . 
All white dwarfs therefore have been traditionally found to belong to one of two 
categories: those with a hydrogen-rich atmosphere (the DA spectral type) and 
those with a helium- rich atmosphere (the non-DAs). Here we report the discov- 
ery of several white dwarfs with atmospheres primarily composed of carbon, with 
little or no trace of hydrogen or helium. Our analysis shows that the atmospheric 
parameters found for these stars do not fit satisfactorily in any of the currently 
known theories of post-asymptotic giant branch evolution, although these objects 
might be the cooler counterpart of the unique and extensively studied PG1159 
star H1504+65 (refs 4-7). These stars, together with H1504+65, might accord- 
ingly form a new evolutionary sequence that follow the asymptotic giant branch. 

Traces of carbon are typically observed as either neutral carbon lines or molecular C2 
Swan bands (defining the DQ spectral type) in cool helium-rich white dwarfs with effective 
temperatures (T e g) below ~ 13,000 K. The presence of carbon in the atmospheres of these 
objects has been explained successfully by a model in which carbon is dredged-up from the 
underlying carbon/oxygen core by the deep helium convection zone 8 . This model predicts a 
maximum contamination of carbon at an effective temperature of ~10,000 K (corresponding 
approximately to the temperature at which the surface convection zone is maximal) be- 
fore gradually decreasing with lower temperature, in agreement with atmospheric analysis 
determinations 9 . We note that although some of these objects show a very high level of 
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carbon pollution, helium always remains the dominant constituent of the atmosphere (the 
highest carbon abundances by number (N) found are around logN(C)/N(He) —3). 

Atmospheric abundance determination for DQ white dwarfs has been most successful 
for stars on the cool side of the maximum contamination and very few analyses have been 
done for hotter objects located on the ascending side of the curve 10-12 . Recently 13 , several 
new hot DQ white dwarfs have been discovered in the Sloan Digital Sky Survey (SDSS), 
providing a unique opportunity to test the dredge-up theory on the "hot side". However, 
appropriate atmospheric models were not available at that time. Because the atmospheric 
compositions of these white dwarfs were expected to be helium-dominated, effective tem- 
perature estimates (Tff e ) were obtained by comparing the SDSS colours with those of pure 
helium composition models. Although it was expected that such an approximation would 
overestimate the effective temperature by several thousand kelvin, some of these stars were 
so hot (Tff e > 25,000 K) that an overlap with most of the helium-rich star (the DB spec- 
tral type) temperature range seemed unavoidable. It was hypothesized that these stars had 
thinner outer helium envelopes so that dredge-up occurred earlier in the cooling sequence. 
These highly carbon-polluted white dwarfs are expected to be massive, so it was argued that 
they might represent the missing high-mass tail of the DB mass distribution 14 . 

Thus, it is with this scientific rationale in mind that we proceeded with the calculation 
of the appropriate atmospheric models for these objects. Because the continuum opacity 
of heavy elements might not be negligible for these objects, these new models have been 
updated with the latest C and O photoionization cross-sections from the Opacity Project 15 . 
Although the analysis of the coolest objects (T c s < 15,000 K) is straightforward (results will 
be presented elsewhere; manuscript in preparation), we quickly realized that no combination 
of carbon and helium could successfully reproduce the observed features (mostly Cll lines) 
in the optical spectrum of the hottest ones by assuming a helium-dominated atmosphere. 
Indeed, such models predict the presence of a strong Hel A=4,471 line that is not observed 
spectroscopically in our sample of hot DQ stars. We thus concluded that a good fit to both 
the spectra and the energy distribution was possible only by considering atmospheres made 
primarily of carbon, with little or no trace of hydrogen and helium. 

Figure 1 shows our fits to the optical spectrum and photometric energy distribution of 
one such star using a grid of pure carbon atmosphere models, as well as models containing 
traces of hydrogen and helium. The fitting method is similar to that used for cooler DQ and 
needs not to be repeated here 9 . This particular case shows a small trace of hydrogen which 
we believe is more likely to be the result of accretion from the interstellar medium than of 
primordial origin. We found eight more similar objects in the SDSS white dwarf catalogue 16 , 
all of which were found to have a carbon-dominated atmosphere and a temperature between 
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18,000 and 23,000 K, although for some faint stars, detailed analysis is more uncertain 
because of the poor signal-to-noise ratio of the observations. Only crude upper limits for 
hydrogen and helium (of the order of logN(C)/N(H) > 1.5 and logN(C)/N(He) > 1.5) can be 
obtained from the other spectra. Our detailed model atmosphere fits provide fair estimates 
of the effective temperature and chemical composition, but higher- quality data (which we 
hope to obtain in December 2007, weather permitting) are essential for a precise measure of 
the surface gravity (and thus mass) of these objects. 

These new stars are too hot to be explained by the standard convective dredge-up 
scenario and so we sought another explanation. Natural progenitor candidates to consider 
are the hot PG1159 stars. The latter are probably the result of a late He-shell flash 17 (the 
so-called "born- again" scenario) at the end of the post-asymptotic-giant-branch phase, which 
almost completely eliminates the remaining hydrogen and mixes material from the interior 
with the helium envelope. As a result, these objects re-enter the white dwarf cooling track, 
but this time they have a surface composition that is a mixture of helium, carbon, oxygen, 
and little or no hydrogen (T e ff between ~75,000 and 200,000 K; typical abundances 18 , in 
mass fraction, of He, C, O and Ne are 33%, 50%, 15% and 2%). However, as these stars cool 
down, gravitational diffusion rapidly separates the heavier elements from the helium that 
tends to float to the surface and by the time they have reached 25,000 K, helium completely 
dominates the surface composition. 

Models exploring the evolution of stars 9-11 M Q have produced white dwarfs with O- 
Ne-Mg cores and CO envelopes and possibly little He and H; such models may explain our 
peculiar stars 19-21 . However, it is unclear what quantity of He and H should remain for 
such stars. Furthermore, the resulting white dwarfs have high surface gravities (log g ~ 9, 
in centimetre-gram-second units) and would thus show extremely broad carbon lines that 
are incompatible with the observed widths of the lines (Fig. 1 shows that log g 8 is 
more likely, although analysis of a better spectrum is needed for more precision). The most 
likely scenario to explain the very existence of hot DQ white dwarfs with carbon-dominated 
atmospheres is that these stars are the progenies of objects such as H1504+65. The latter 
is a unique object among the known PG1159 stars so far. It is the hottest specimen of 
its class at T cS 200,000 K, and its atmospheric composition very unusual, with a mass 
fraction of ~ 50% C and ~ 50% O plus small traces of heavier elements, but no detectable 
helium or hydrogen. It is currently believed that HI 504+65 is essentially a bare stellar 
nucleus produced by a particularly violent post-asymptotic-giant-branch very late thermal 
pulse that has destroyed in large part the remaining stellar envelope containing helium and 
hydrogen. The best available simulations of such an event, those of ref. 17, suggest indeed 
that, following the last He flash, the full envelope becomes convective and extends to deep 
enough layers for H and He to be completely consumed, except that very small traces of He 
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may survive. 

The uniformity of the chemical composition of the envelope of H1504+65 is believed to 
be maintained for a while by a residual wind but, with time, diffusion becomes dominant and 
He, C and must separate under the influence of gravitational settling. The idea is that 
the wind slowly dies out with cooling (decreasing luminosity) and that the braking effect it 
provides on element separation becomes less and less efficient. Hence, the small residual trace 
of helium believed to exist in the envelope of H 1504+65 eventually diffuses upward to form 
a thin layer above a C-enriched and O-depleted mantle. The total mass of an atmosphere 
is tiny (~ 10~ 14 to 10 -15 M*, where M* is the mass of the star), so there is ultimately 
enough accumulated He to form a full atmosphere and the descendant of HI 504+65 would 
now be "disguised" as a He-atmosphere white dwarf after the PG1159 evolutionary phase. 
With further cooling, a convection zone develops in the C-enriched mantle owing to the 
recombination of that element while the overlaying He layer remains radiative. At some T e g, 
whose exact value is not yet known due to a current lack of proper models, it is conjectured 
that the subphotospheric C convection zone becomes active enough to be able to dilute from 
below the overlaying He layer. At that point the star would undergo a dramatic spectral 
change, transforming itself from a He-dominated atmosphere white dwarf to a star with a 
C-dominated atmosphere, because the mass in the C convection zone is orders of magnitude 
larger than the mass of the He layer. 

Hence, the former PG1159 star H1504+65, showing initially a mixed C and O atmo- 
sphere, would now show a C-dominated atmosphere after an intermediate phase in which it 
would have been observed as a He-atmosphere white dwarf. We note, however, that helium 
must ultimately reappear at the surface in enough quantity to form a helium-rich atmosphere 
because no carbon-rich object has ever been discovered at lower temperatures (T eff < 15,000 
K). The remaining helium probably floats again at the surface when the star cools down, 
perhaps turning the carbon-rich objects into DQ stars belonging to the second sequence of 
high carbon abundances 11 ' 22 . We are at present making evolutionary calculations to test the 
above scenario. 

To conclude, we roughly estimate that the space density of carbon-rich white dwarfs, 
assuming a 50% completeness for the targeting in SDSS 16 , is between ~ 2.1 x 10~ 6 pc~ 3 
and ~ 7.0 x 10~ 8 pc~ 3 depending on the preferred value of log g (and thus radius) and T c s 
for the stars. For comparison, the local space density of white dwarfs 23 is estimated to be 
~ 5.0 x 10~ 3 pc~ 3 . White dwarfs with carbon dominated atmospheres are thus intrinsically 
rare and a relatively large volume of space (as was done by the Sloan Digital Sky Survey) 
had to be surveyed before many of them could be found. 
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Fig. 1. — Fit to the optical spectra and energy distribution for a carbon-rich white 
dwarf, a, The thick red line represents our log g = 8 best fit (the parameters are indicated 
in the panel) while the blue line represents a solution for logg = 9. Hydrogen abundance is 
determined by fitting H/3 (A=4,86lA) while helium abundance is constrained by the absence 
of the Hell A =4,47lA line, b, Photometric measurements in the u, g, r, i and z bands are 
represented by error bars, while the average model fluxes from the same best log g = 8 model 
are shown by filled circles. Both a and b have the same y axis. 
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